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Abstract 

The use of Pade approximants for the description of QCD matrix elements is discussed in this talk. We will see how 
they prove to be an extremely useful tool, specially in the case of resonant amplitudes. It will allow the inclusion 
of high-energy Euclidian data to improve the determination of low-energy properties, such as the quadratic vector 
radius. This does not mean that the rational approximations can be arbitrarily employed for the extraction of any 
desired hadronic parameter. A discussion about the validity, limitations and possible issues of the Pade analysis is 
carried on along the paper. Finally, based on the de Montessus de Ballore's theorem, a theoretically safe new procedure 
for the extraction of the pole mass and width of resonances is proposed here and illustrated with the example of the 
p(770). 
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1. Introduction 

Quantum Chromodynamics (QCD) has been proven to be the right theory to describe the strong dynamics interac- 
tions. However, in the non-perturbative regime one finds resonant structures in the cross section which can be related 
with the presence of short-lived intermediate hadronic states, usually referred as resonances. From a mathematical 
point of view, these states appeared as complex poles of the amplitude in the transferred energy at higher complex 
Riemann sheets, other than the physical one. Typically, the amplitude grows abruptly when the energy gets close to 
the resonance pole, remaining nevertheless finite all the way as the singularity is located off the real axis, at a complex 
value of the energy. 

One of the clearest examples of this kind of resonant amplitudes is the 7T7r-vector form-factor (VFF). There, the 
spectral function is dominated by a very pronounced p(770) meson peak and essentially no other big effect is observed. 
Through this qualitative knowledge of the spectral function in the Minkowsky region (q 2 > 0) and the experimental 
data from the Euclidean region (q 2 < 0) one is able to extract properties of the VFF at the origin 1 1 ] . Its first and second 
derivatives were determined at q 2 = (respectively related to the quadratic vector radius ( r 2 )y and the curvature c" v ) 
by means of Pade approximants (PA) centered at the origin | lj]. However, this precise procedure does not allow us to 
make predictions for properties of the amplitude above threshold, such as the p-meson pole position. 

Indeed, the potential danger of using rational approximants for the extraction of resonance pole positions is shown 
with the help of a model. Under some limits the PAs result equivalent to some unitarization procedures, such as the 
inverse amplitude method (IAM) J^EL reason why these unitarizations have been sometimes loosely called "Pades". 
We performed a perturbative computation in the linear sigma model (LSM) and found that the PA around the origin 
led to improper determinations of the meson mass and width [4]. However, we will show here that by an adequate 
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reinterpetation of the PA it is possible to converge (though slowly) to the actual low-energy constants (LEC) of the 
model. Nonetheless, the perturbative calculation only makes sense in the LSM when the ct-jui interaction is weak. It 
has been argued that the unitarization procedures such as IAM are expected to work for cases of strongly interacting 
mesons, such as the physical cr(600), and it is not intended for weakly interacting theories [5]. This issue is still 
unclear and will require of further clarifications. 

Still, the Pade approximants allow us to produce a model independent determination of the resonance poles if they 
are adequately employed. To do this in a theoretically safe way, we need to center our Pade above the branch-cut 
singularity (beyond the first production threshold), not bellow (at the origin). This modification also makes possible 
the direct use of Minkowsky data, being now the Euclidean ones discarded. 

The importance of these hadronic pole parameters is that usually one relates an observable with the corresponding 
renormalized couplings but in the resonance case in QCD there is still plenty of debate about which is the right 
lagrangian formulation. Alternatively, the resonance pole positions in the complex energy plane are universal for 
all the processes with those same quantum numbers. They do not depend on a particular lagrangian realization. 
Nonetheless, in many cases extracting these hadronic properties brings along much model dependence as it is not 
clear how to extrapolate from the data on the real energy axis into the complex plane. This is highly non-trivial as one 
can see, for instance, observing the broad spreading of predictions for the cr(600) meson pole (I = J = channel in 
7T7r-scattering) [6]. 

The quadratic radius and curvature of the 7T7T-VFF will be extracted in Sec. 2 with the help of Pades centered at 
q 2 = 0. A discussion on Pade unitarizations is provided in Sec. 3. Finally, a new kind of PA is proposed in Sec. 4 for 
the study of resonant amplitudes, centered at energies q 2 over the first production threshold. 

2. Pade approximants and the space-like VFF 

Our goal in this section is the description of the 7T7T-VFF F(Q 2 ) in the space-like region: 

<jtV)I I wfu - ~ 1-fd - ~ jy* k»> = (P + p'f m 2 ) , (i) 

where Q 2 = -q 2 = -(p' - p) 2 , such that Q 2 > corresponds to space-like data. Since the spectral function for the 
corresponding dispersive integral for F(Q 2 ) starts at twice the pion mass, the form factor can be approximated by a 
Taylor expansion in powers of the momentum for \ Q 2 \ < (2m n ) 2 . 

We want to construct an approximation that can be systematically improved upon. However, it will not be our aim 
to extract time-like properties from this analysis, like, for instance, vector meson mass predictions. It is neither our 
intention to describe the amplitude on the physical absorptive cut, above the jtjt threshold. Finally, it is convenient 
to remark that the results [ 1] presented here do not refer to any large-Afc approximation but to the physical Nc = 3 
quantities. 

2.1. The method: Pade approximants 

Consider an analytical function F(z) at a point, e.g., z — 0. A Pade approximant P%(z) = 6iv(0/^Af(z) is defined 
by the ratio of two polynomials of degrees N and M which agree with the original function F(t) up to the derivative 
of order N + M at z = 0: 

P N M (z) - F(z) = ^(z™ 1 ). (2) 

One may wonder what is new here with respect to a Taylor series of the form F (z) — «o + o-\Z + atz 2 + ... The 
difference is that the polynomials are unable to go beyond the singular points. They are not able to describe them, 
setting those singularities the maximum size of the convergence disk (centered at the analytical point z - 0). On the 
other hand, the poles of the PA tend to mimic the singular structure of the original function F(z). For instance, if one 
studies the P^,(z) approximant of In (1 + z) all the poles are generated at the branch cut -oo < z < -1, which gets more 
and more densely populated as N — > oo. 

Thus, in many cases, the PAs are found to work far beyond the analytical disk of convergence of the Taylor series. 
This allowed us to use space-like data to improve our description of the VFF at the analytical point Q 2 = 0. However, 
in general, the Pades converge over a compact region and, although this region may get larger and larger as one 
increases the order of the Pade, we did not use in this analysis information from the VFF at Q 2 = oo. 
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It is convenient to remark that in the case of non-meromorphic functions, such as the physical amplitudes with 
logarithmic branch-cuts, the poles of the PA centered at the origin do not correspond to resonances, but rather to 
bumps and other structures in the spectral function over the logarithmic branch cut, which the PA tends to mimic. 

In our phenomenological analysis of the VFF, we will use as inputs all the available data in the Euclidean region, 
which range from Q 2 = 0.01 GeV 2 up to 10 GeV 2 . We will also make use of the qualitative knowledge we have on 
7T7T-VFF spectral function p(s), essentially provided by the p(770) peak. This suggests the use of the sequence for 
the description of this particular observable. The final aim of the analysis will be the extraction of the first and second 
derivatives at q 2 = 0, i.e., the quadratic vector radius ( r 2 ) n v and the curvature c n v . From this perspective, the vector 
meson dominance expression F(Q 2 ) — (1 + Q 2 /M 2 )~ l is just a Pade, the first term of a sequence. 

2.2. Theoretical uncertainties: playing with a phenomenological model 

In order to illustrate the usefulness of the PAs as fitting functions, we will first use a phenomenological model as a 
theoretical laboratory to check our method. The model will also give us an idea about the size of possible systematic 
uncertainties. 

We will consider a VFF phase-shift with the right threshold behavior and roughly the physical values of the rho 
mass and width. The form-factor is recovered through a once-subtracted Omnes relation, 



F(Q 2 ) = exp 



T J 4, 



6(f) 



dt 

4ml t(t + Q 2 ) 



(3) 



where 8{t) plays the role of the vector form factor phase-shift [7, 8, 9] and is given by 



5(f) = tan 



MSJt) 



Mj-t 



(4) 



with the f-dependent width given by 



r P (0 = r 



M 2 
\' v1 p > 



e(t-4ml) 



(5) 



and cr(t) - J 1 - 4m 2 / 1 . The input parameters are chosen to be close to their physical values: Tq — 0.15 GeV, 
Mp = 0.6 GeV 2 , 4m 2 = 0.1 GeV 2 . This model is actually quite realistic and it has been used for the extraction of the 
physical mass and width of the p(770) meson from time-like experimental data 10, [H @] • 

We generate now an emulation of the experimental data from our theoretical model. In order to recreate the 
situation of the experimental data lfloll - lfl5ll with the model, we have generated fifty "data" points in the region 0.01 < 
Q 2 < 0.25, thirty data points in the interval 0.25 < Q 2 < 3, and seven points for 3 < Q 2 < 10 (all these momenta in 
units of GeV 2 ). These points are taken with vanishing error bars since our purpose here is to estimate the systematic 
error derived purely from our approximate description of the form factor. 

These generated data is then fitted through P\ Pade approximants, 



P\(Q 2 ) = l + J]a k (-Q 2 ) k +(-Q 2 ) L 



a L 



(6) 



where the vector current conservation condition Pf(0) = 1, i.e., ao = 1, has been imposed. At low energies this 



produces for the Taylor coefficients aj the prediction 

F(Q 2 ) = 1 - a, Q 2 + a 2 Q 4 - a 3 Q 6 + ... 



(1) 



This leads to a series of predictions for the low-energy parameters, which are compared to their (known) exact 
values in TableQ] The last PA we have fitted to these data is P b { . Notice that the pole position of the Pades differs from 
the mass parameter of the model M 2 -and this from the pole mass-. This example makes explicit how one can get a 
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PI 


P\ 


P\ 


P\ 


P\ 




P\ 


F(g 2 )(exact) 


ai(GeV- z ) 


1.549 


1.615 


1.639 


1.651 


1.660 


1.665 


1.670 


1.685 


a 2 (GeV- 4 ) 


2.399 


2.679 


2.809 


2.892 


2.967 


3.020 


3.074 


3.331 


a 3 (GeV- 6 ) 


3.717 


4.444 


4.823 


5.097 


5.368 


5.579 


5.817 


7.898 


s p (GeV 2 ) 


0.646 


0.603 


0.582 


0.567 


0.552 


0.540 


0.526 


M 2 p = 0.6 



Table 1: Results of the various fits to the form factor F(Q 2 ) in the model, Eq. The exact values for the coefficients a, in Eq. (Q 
are given on the last column. The last row shows the predictions for the corresponding pole for each Pade (s p ), to be compared to 
the mass parameter Mr p = 0.6 GeV 2 in the model. 



ai (GeV _i ) 



a 2 (GeV ) 



Figure! : 

fits fion- 



ue of the a\ and ai Taylor coefficients for the P\ sequence of Pade Approximants obtained from experimental data 



rather precise value for the Taylor coefficients in Eq. © without an accurate knowledge of the spectral function (i.e., 
of the time-like region). 

Based on the previous results, we will take the values in Table[T]as a rough estimate of the systematic uncertainties 
when fitting to the experimental data in the following sections. Since, as we will see, the best fit to the experimental 
data comes from the Pade P\ , we will take the error in Table[T]from this Pade as a reasonable estimate and, respectively, 
add to the final error an extra systematic uncertainty of 1 .5% and 10% for a\ and a 2 . 

2.3. Experimental pion vector form factor 

The prominent role of the rho meson contribution motivates the use of the P\ Pade sequence as the central tool for 
the study of this amplitude, later complemented by other types of Pades. The fit of P\ to the space-like data points in 
Refs. [ 10]— [ 15] determines the coefficients a^ that best interpolate them. Fig.[T]shows the evolution of the fit results 
for the Taylor coefficients a\ and a 2 for the P L X PA from L — up to L — 4. As one can see, after a few Pades these 
coefficients become stable. For the data in Refs. H 1 Ofl — HI 5fl ^ this happened at L = 4. The P 4 Pade Approximant 
provides our best fit and, upon expansion around Q 2 = 0, this yields 

a x = 1.92 ±0.03 GeV" 2 , a 2 = 3.49 ± 0.26 GeV" 4 ; (8) 

witha^ 2 /dof = 117/90 H. 

Eq. © shows that the pole of each P\ PA is determined by the ratio s p = a^ja^x. This ratio is shown in Fig. [2] 
together with a gray band given by M 2 + M p T p for comparison. From this figure one can see that the position of the 
pole of the PA is close to the physical mass M 2 [6], although it cannot be identified with it, as we already saw in the 
model of the previous subsection. The P\ pole s p is always real and lies at the "bump" of the spectral faction, in the 
range M p + M p T p . The Pade tends to reproduce the p peak line-shape but, obviously, no complex resonance pole can 
be recovered from a Pade at higher Riemann sheets. 

As one can see in Fig. [3] the sequence converges to the physical form-factor in the data region but, eventually, 
it diverges like (Q 2 ) i_1 . These PAs only converge on a compact region of the complex plane and are unable to recover 
the 1 IQ 2 asymptotic behaviour prescribed by QCD at short distances 01611 . Nonetheless, the important fact is that the 
Pades allow the use of data not only by the origin but even from energies as large as Q 2 ~ 10 GeV 2 , something that a 
normal Taylor expansion at the origin does not permit. 

This PA analysis was complemented with the results from other types of Pades: two-pole Pade approximants 
PA\\ one, two and three-pole Pade-types (PT) with the poles fixed beforehand; two-pole partial Pades PP\ X , with one 
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s p (GeV 2 ) 

0.7 ; 



Figure 2: Position s p of the pole for the different Pj. The range with the physical values M 2 ± M p T p is shown (gray band) for 
comparison. 



Q'F(Q 2 ) (GeV ) 



" g 2 (GeV) 

2 4 6 8 10 12 

Figure 3: The sequence of P\ PAs is compared to the available space-like data [ 10] -1 15]: P\ (brown dashed), P\ (green thick- 
dashed), P 2 (orange dot-dashed), Pj (blue long-dashed), P\ (red solid). 





PA 2 


PTf (P) 


pt% (p,p0 


PT% ip,p") 


PT£ (p,p',p") 




ai(GeV- 2 ) 
a 2 (GeV- 4 ) 


1.924 ±0.029 
3.50 ±0.14 


1.90 ±0.03 
3.28 ± 0.09 


1.902 ±0.024 
3.29 ±0.07 


1.899 ±0.023 
3.27 ± 0.06 


1.904 ±0.023 
3.29 ±0.09 


1.902 ±0.029 
3.28 ±0.09 



Table 2: Different results for two-pole Pade approximants, Pade-types and partial Pades, where we used the p, p', and p" masses. 



pole fixed beforehand and the other determined by the low-energy coefficients a^. The different results are gathered 
in Table. |2 

Combining all the previous rational approximants results in the average given by |l]] 

a x = 1.907 ± 0.010 s ta t ± 0.03 syst GeV" 2 , a 2 = 3.30 ± 0.03 stat ± 0.33 syst GeV" 4 . (9) 

The first error comes from combining the results from the different fits by means of a weighted average. On top of 
that, we have added what we believe to be a conservative estimate of the theoretical (i.e. systematic) error based on 
the analysis of the VFF model in the previous subsection. We expect the latter to give an estimate for the systematic 
uncertainty due to the approximation of the physical form factor with rational functions. For comparison with previous 





(r% (fm 2 ) 


a 2 (GeV- 4 ) 


This work 


0.445 ± 0.002 stat ± 0.007 syst 


3.30 ± 0.03 stat ± 0.33 syst 


PDG [6] 


0.452 ±0.011 




CGLIT7, 18J 


0.435 ± 0.005 




TY [19] 


0.432 ±0.001 


3.84 ±0.02 


BCT [20] 


0.437 ±0.016 


3.85 ±0.60 


PP [9] 


0.430 ±0.012 


3.79 ±0.04 


Lattice [2J_] 


0.418 ±0.031 





Ta b le. 3 : O ur re sults for the quadratic vector radius (r 2 Y v and second derivative «2 are compared to other determinations I6l.l9l.fi7l 
|ialia|20|,|2l|]. Our first error is statistical. The second one is systematic, based on a previous analysis of a VFF model. 
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analyses, we also provide in Table[3]the value of the quadratic vector radius, which is given by {r 2 ) n v — 6 a\ . 

In summary, we used rational approximants as a tool for fitting the pion vector form factor in the Euclidian range. 
Since these approximants are capable of describing the region of large momentum, we think they are better suited 
than polynomials for the description of the currently available space-like data. As our results in Table 2 show, the 
errors achieved with these approximants are competitive with previous analyses existing in the literature, based on 
more elaborated techniques. 



3. A critical look on Pade unitarizations 



Exact unitarity is an important additional piece of information that eventually needs to be included in the de- 
scription of the scattering processes. A commonly employed unitarization procedure is the inverse amplitude method 
(IAM) 10,13, 

l22l |23L 12411 . Although it can be formulated in a more elaborated way through dispersion relations, the 
method relies in the unitarity relation in the elastic region, which in the massless case has the simplified form, 



lmt(s) = \t(s)\ 2 



lm[t( S y l ] = -1 



(10) 



This fixes completely the imaginary part of the inverse partial-wave amplitude t(s)~ l on the elastic part of the right- 
hand cut (s > 0). All that remains in f(s) _I = Re[t(s)~ l ] + ImfX,?) -1 ] is to determine the real part of t(s) , which is 
fixed in the IAM through a low-energy matching to ^PT: 



K s )xPT - f (2) + f (4) + *(6) + 



(11) 



where the t^) are the contributions corresponding to Q(p k ) in^PT. Hence, depending on the order of the matching one 
obtains a sequence of unitarized amplitudes QE, 22, 23, 24]: 



0( P 4 ) 



tlAM - 



1 - f(4)A(2) 



0( P 6 ) 



tlAM - 



1 - f (4)/ f (2) - f (6)/f(2) + (f(4)A(2)) 



2 ' 



(12) 



~ 1 GeV . 
/o(980)Si 



Identical results are obtained if one recovers the partial wave amplitude through a dispersion relation and matches 
^ PT on the left-hand cut. 

In the tree-level limit, the IAM expressions ( fT2l become a series of Pade approximants of the form P\, Pi, etc, 
reason why these unitarizations are sometimes called "Pades". In any case, it is in this limit that one can use all the 
powerful technology of the mathematical theory of Pade approximants. However, it has been argued that the IAM 
should be only applicable for the description of broad resonances, such as the sigma meson, and not narrow states . 
The IAM has been found to described the data reasonably well even, in some cases, up to energies as high as 
|3fl . It has been able to generate all the resonances below 1 GeV: p, K*,o~ (or /o(600)), k, ao(980) and 
23TI24I1 . and the expected Nc behaviour of the p and K* poles as qq states has been recovered I123ll . 
However, one may wonder what information is lost when Re[f(i) _1 ] is fixed at low energies with^PT. Likewise, in 
general the unitarized amplitudes violate crossing, as only it only resums a particular set of diagrams in the s-channel, 
not in the crossed ones. Nevertheless, even though the IAM determinations have produces good numerical results, 
this is not supported by the theory of Pade approximants. Thus, we consider that the theoretical reason why the IAM 
has been so successful at the phenomenological level still needs and deserves further clarifications. 

3.1. A counterexample: The cr in the Linear Sigma Model 

The application of the IAM to the linear sigma model (LSM) was proposed in Ref. [4] as a counterexample, to 
show how the prediction from unitarizing the low-energy LSM led to very different pole mass and width than those 
actually in the LSM. 

At tree-level, the sigma mass is found to be M 2 = 2p 2 , and the width is zero. At next-to-leading order, the sigma 
pole gets shifted due to the quartic potential, i.e., = 2/j 2 + 0(g), and the width becomes different from zero. For 
simplicity, the massless pion limit is assumed. 

In order to determine the scalar meson mass and width up to 0(g), we compute the one-loop sigma correlator [25], 



iA(sy 



Mt 



1 + 



3g 
16tt 2 



13 

T 



In 



Mi 



+ 3p(s) In 



P(s) + 1 
Pis) - 1 



+ 0(g 2 ) 



(13) 



where p(s) = Jl — 4M*/s and the term -13/3 is determined by the renormalization scheme chosen by Ref. 12511 . 
which sets the relation 2gF 2 = M\ at the one-loop order, with F the pion decay constant and M a the renormalized 
mass parameter. Now it is possible to extract the pole s p of the propagator up to the considered order in perturbation 
theory. If one approaches the branch cut from the upper part of the complex s-plane, the pole in the second Riemann 
sheet is located at s p = (M p - iF p /2) 2 , with the pole mass and width, 



p 
M 2 



1 + 



167T 2 



13 

y 



■n^+0(g 2 ) 



M p r„ 

Ml 



The 7T7r-scattering is determined by the n + n 



T(s, t, u) 
T(s, t, u) 



1=0 



i=i 



1=2 



T(s, t, u) 

and the partial wave projection provided by 



3s 



+ 0(g 2 ). 



I6n 

amplitude A(s, t, u). This defines the isospin amplitudes 



(14) 



3A(s, t, u) + A(t, s, u) + A(u, t, s) , 
A(t, s, u) - A(u, t, s) , 
A(t, s, u) +A(u, t, s) , 



1 f 1 

f 7 ,0) = d cos 6 Pj(cos 8) T(s,t,uf , 

J 64ttJ _j 



(15) 



(16) 



where 8 is the scattering angle in the nn center-of-mass rest frame. 

We now consider the LSM at low energies, which reproduces the structure prescribed by ^PT. Hence, for the first 
partial waves tj(s), with I J = 00, 1 1, 20, it produces the 0(p 2 ) amplitudes, 



and at 0(p A ), 



2 ' 



\6nF 



<o(*)(?) x 



96ttF 



2 ' 



t Q (s)(2) = ~ 



32nF z 



lis 
6M 2 



264n z 



-s s 193 

18 In— l+71n-- + — 



Mi 



Mi 



+ 0(g 2 ) 



(17) 



—s 

M 2 

-2s 
3M! 



1 + 



487r z 



Ml Ml 3 1 



24t!- 2 \4 Ml 



-s 11, s 163 
In — - + — In — - + 



Mi 



24 



+ 0(g 2 ) 



(18) 



The Inverse Amplitude Method (IAM) provides an amplitude that is unitary not only at the perturbative level but 



exactly. At 0(p ), one has the unitarized amplitude, t u 



fa) 



f(4) 



which has its poles s p at t(2)(s p ) = t^)(s p ): 



t{2) 



I J =00 

IJ=11 

IJ=20 



s„ = 



6 2 
— M 
11 ° 



-M 

3 



1 + 



2647T Z 



193 6 i 

— +25 In— -mn\ + 0{g 2 ) 



1 + 



= --Mi 



48^- 
1 + 



26 



+ in\ + 0(g 2 ) 



8 /163 



Hot 



24.M24 +51n 2 + 4 



+ 0(g 2 ) 



(19) 
(20) 
(21) 



These are the poles that appear in the unphysical Riemann sheet as one approaches from upper half of the first Riemann 
sheet. There is also a conjugate pole at s* p if one approaches the real s-axis from below. 
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The first thing to be noticed is that poles appear in the IJ — 1 1 and 20 channels even for small values of g, contrary 
to what one expects in the LSM, where no meson with these quantum numbers exists. Furthermore, these "states" are 
not resonances, as they are located on the left-hand side of the complex s-plane, out of the physical Riemann sheet, 
and carrying a negative squared mass. 

One can easily see the important disagreement between the actual LSM value for the pole mass and width given 
in Eq. ( TBI and the posterior IAM "prediction", 



1 M 2 ^ 
M 2 

1V1 IT ! 



_6_ 
IT 



1 + 



* /50 



In- 



167r 2 \33 11 



386 
~99~ 



+ 0(g 2 ) 



Ml 



121 I6n S ' 



(22) 



this is, (M 2 ,) IAM - 50% (M 2 p ) LSM and {M p T p ) IAM - 20% (M p T p ) L 



3.2. Higher order Pad Approximants for the LSM 

In the tree-level limit (for instance, at large Nc), the IAM amplitudes become Pade approximants centered 
at s = 0: 



[<V)1 



t(2) 



f (2) + ?(4) + ■■ 



10(P 6 )] 



1 - ?(4)/?(2) - f (6)/f(2) + (?(4)/?(2)) 



1 - ?(4)/f(2) 

t(2) + f (4) + ?(6) + — 



Pl(s) 



The 7T7r-scattering is given in the LSM case by 



s M 2 

A(s,t,u) = —? — r 

F 2 Ml - s 



«o(*) 



Ml 



32ttF 2 



3M± 2Ml ( s 

-5+ ' T + -In 1 + — r 

Ml-s s \ Ml 



(23) 



which at low energies become, 



A(s,t,u) = — 



1 + 



Ml 



Ml 



t° (s) = 



16jtF 2 



1 + 



Us 
6M 2 



15/ 
l2Mt 



(24) 



The first Pade approximant, P\, gives for the <x pole the prediction 



= tjM 2 . But how does the series P\, P\, 

We 



P\... evolve? One could even wonder, for instance, about the behaviour of other sequences such as P\, P\, Py 
have shown in Fig. [4] the P l M sequence up to such a high order as P l 6y One can see the extremely slow converges, 
with the pole prediction still a 30% off for P l 6l . Furthermore, this kind of P l M Pade approximants has only one zero (at 
5 = 0) and places the M poles in a circle centered at s = 0, producing an analytical structure completely different to 
that in the actual LSM. 

On the other hand, a quick convergence is found for the diagonal Pade sequence P^: P\ reproduces the sigma 
pole a 40% off but P\ disagrees by less than 1%, P\ by less than 0.1%, etc. Likewise, Fig.[5]b shows how the P^ 
sequence, besides providing the isolated pole of the sigma, tends to reproduce the left-hand cut as increases. The 
poles of P^J are plotted there. Although a PA is a rational function without cuts, these are mimicked by placing poles 
where the cuts should lie. The P 2 ^ has one isolated pole near Ml (with an accuracy of 10~ 30 ) and nineteen poles over 
the real axis at s p < -Ml, i.e. on the left-hand cut of the original function. As is increased, the number of poles 
lying on the branch cut increases too. A last remarkable feature is that the P^ approximants obey exact unitarity, as it 



happened with the IAM sequence P N . 

However, phenomenologically, the IAM has been usually employed for the study of the experimental data in the 
resonance region. These have been used to fit its parameters (e.g., the (i in Ref. yfl) and to extract the resonance 
pole positions. These parameters also provide a prediction for the chiral low-energy constants. However, although in 
general the convergence of these parameters to the LECs is unclear, in the case of the sequence of one-pole Pades P^ 
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Figure 4: (a) Position of the nearest pole to M\ for the first PAs of the form P [ N with N odd (for even N all the poles are complex), 
(b) Poles of the PL in the complex plane. 
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Figure 5: (a) Location of the closest pole to for the first PjJ Pade approximants. (b) Poles of Pf' y 



centered at s ~ M~ , the convergence is ensured in the disk centered at that point and with maximal radius limited by 
the position of the essential singularity of the left-hand cut ats = —M\ 1261, 12711 . Hence, the predictions of the P N X Pade 
approximant at the middle point 5 = converge to the actual ^PT low-energy couplings: the P\ centered at s = 
recovers the 0(p 4 ) coefficient with a 50% error, 31% for P 2 V 17% for P\, etc. From this perspective, maybe one could 
explain the phenomenological success of the IAM (which shares the first term, P\), although the P^ sequence might 
be more adequate for low-energy predictions than the usual IAM pattern PL 



4. Model independent determination of the resonance poles 

As we have seen in previous sections, the Pades have been used to extract low-energy parameters from Euclidean 
data. In a missleading way, some unitarizations procedures have been called Pades, although their justification from 
the point of view of Pade theory is unclear. Though these unitarizations have been proven more or less successful in 
their predictions of the resonance poles and low-energy coefficients, they are not supported by the mathematical theory 
of Pades. Alternatively, in this section we show how it is possible to construct a different type of Pade approximants 
which allows us to use Minkowskian data and to extract the resonance pole position through a theoretically safe 
procedure supported by mathematical theorems. 

To illustrate the procedure, let us start by the simplest possible case. If one has a function F(s) analytical in a disk 
Bg(so) then the Taylor series 5jy(s) = Yjk=o a k( s ~ s o)* converges to F(s) in Bg(so) for N — > oo, with the derivatives 
at = F^ k '(so)/kl. Experimentally, one usually does not have the derivatives at some point, s = so, but a series of 
experimental points Fj at different s } ; from which one extracts the function and its derivatives through polynomial fits 
with S n(s) at higher and higher order N. 

But, what happens if there is a single pole at s = s p within the disk and F(s) is analytical everywhere else in 
Bs(so)7 In that case, the Taylor series does not converge any more. Nonetheless, the needed modification is not 
really big. In its simplest version with just one pole, the de Montessus de Ballore's theorem |26[ 27] states that the 



sequence of one-pole Pade approximants f"/ around sq converges to F(s) in any compact subset of the disk excluding 
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Figure 6: a) Analytical structure of the VFF bellow the first production threshold, b) Structure and analytical extension of the 
1" Riemann sheet between the first and second production thresholds, with a resonance pole s p in the 2" d Riemann sheet in the 
proximity of Sq. 



the pole s p : 

jv-i , _ ,N 
P»(s;s ) = ^> (*-*,)* + aZ ■ ( 25 ) 

k=0 1 - (s - Sq) 

on 

Hence, one finds that the Pade pole x p - so + converges to s p for N — » oo. Experimentally, as referred before, 
one is not provided with the derivatives F(sq), F'(sq)-- but with the values Fj at different Sj, We use then the rational 
functions as fitting functions (in the way done before with the polynomials). As N grows gives an estimate of 
the series of derivatives { F^ k \so) } and, hence, of the pole position s p . 

Usually, the Pades have been constructed around the low-energy point Sq = (with s, typically the total square 
momentum). In matrix elements F{s) without left-hand cut, the amplitude is analytical from s = -oo up to the first 
production threshold s,h and within the disk B„ th (0) (see Fig.|6]a). For instance, one has s,h = 4m% in the nn vector 
form-factor case. Experimentally, one may then have Euclidean data F r (s) at s < and use them to extract the 
derivatives of the VFF at s — OJJ] . Likewise, one may have Minkowskian data F LV (s + i0 + ) from s > s t h which, 
strictly, cannot be used by Pades centered at so = due to the essential singularity at s — s,/,. 

Alternatively one can use in a theoretically safe way Pades centered at so + '0 + over the brunch cut between the 
first and second production thresholds, with s t h < so < ~s t h ( see Fig- 0b). In the 7T7T-VFF this would correspond to the 
range between s t h = 4/«^ and T,/, = 4w|, if multipion channels are neglected. One has then an analytical extension of 
our amplitude F(s) in the I st Riemann sheet at s + i0 + into the 2 nd Riemann sheet. Notice that the function F (sq + i0 + ) 
and the derivatives in the parameters are now complex numbers. 

In the case of resonant amplitudes, a single pole appears in the second Reimann sheet in the neighbourhood of 
the real s axis, which can be related to the existence of a hadronic state -resonance- with the quantum numbers of 



that channel. One can use then the de Montessus de Ballore's theorem B26l 12711 for the description of the data in 
the maximal disk shown in Fig. |6]b. If the resonance pole lies within the disk the Pade approximants allow its 
determination in a model independent way. 

4.1. Testing the method through various models 

We consider a series of p-like models of the nn vector form-factor, with a single pole in the second Riemann sheet 

at s p - (p.77 - 5O.I5) GeV 2 and a logarithmic branch cut (starting at s — for sake of simplicity). The considered 
models were 
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Figure 7: a), b) and c): respectively, the convergence of the P}' poles to the "physical" value in the models A, B and C. The distance 
from the fitted pole to the "physical" one is described by dist = [(M 8 ' - M p ) 2 + (T Rt - T^) 2 ] 2 . 



s In — j 

Model C) F(s) = M_ (26) 

(M - iT/2) 2 - s 

with M and F conveniently tuned in each case to produce the pole at s = s p . 

In order to simulate the physical situation we take the model (A for instance) and generate a series of "data" points 
with zero error, which would represent an ideal experimental situation where all the uncertainty would be theoretical. 
We fit the "data" for the modulus and phase-shift of F(s) and extract the optimal complex parameters for each P^ 
Pade. Notice that this does not mean to fit \F{s)\ (or the VFF phase) with a P 1 ^ Pade. The modulus and phase-shift of 
the data are, respectively, fitted with the modulus and phase-shift of P^. The Pade pole 5°' = (M*' - iF r "/2) 2 is found 
to converge to the "physical" s p = (M p - iT p /2) 2 of the model when — > 00. The approaching of the complex Pade 
pole to its limit value can be observed in Fig. [7] It is clear that the rate of converge depends on the kind of model, 
being faster for model A and slower for B and C. 

4.2. Application to experimental data 

We proceed now to analyze the final compilation of ALEPH 7T7t vector form-factor data for the squared modulus 
\F jm (q 2 )\ 2 [28] and the / = J = 1 nn scattering phase-shift 6„„, identical to the nn vector form-factor phase-shift in the 




760 762 764 766 768 77(1 



M p (MeV) 

Figure 8: 68% CL regions for the rho pole mass and width from the different fits. The smallest (cyan) ellipse provides the 
prediction for N = 3 and the following growing orders in N are given by the ovals with larger and larger size. 
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elastic region Am 2 < q 2 < 4m 2 , (if multipion channels are neglected). This will be the range of application of P 1 ^ Pade 
analysis. For N > 3 the fit^ 2 already lies within the 68% confidence level (CL) and becomes statistically acceptable. 
Their corresponding 68% CL regions for the pole mass and width predictions are shown in Fig. [8] The regions from 
the different fits overlap each other in a compatible way. The allowed ranges become larger and larger as N grows and 
the fit contains more and more free parameters. 

At this point one needs to reach a compromise. On one hand the experimental (fit) errors have an statistical origin 
and increase as one considers higher order Pades P^, with a larger number of parameters. On the other, the systematic 
theoretical (Pade) error decreases as N increases and the Pade converges to the actual VFF In the present work we 
have taken N = 6 as our best estimate as the new parameters of Pades with N > 7 turn out to be all compatible with 
zero, introducing no information with respect to P?. Furthermore, the different models studied before show that in 
any case the theoretical errors for mass and width result smaller than 10 '—10 2 MeV for N > 6, being negligible 
compared to the 0(1 MeV) experimental errors. This yields the determinations 



M n 



763.7 ± 1.2 MeV, 



144 ± 3 MeV , 



(27) 



which is found in reasonable agreement with former determinations obtained from more elaborated procedures and 
with similar size for the uncertainties: 



[Ananthanarayan et al. 112911 ] 
[IAM Q ] 



[Zhou et al. fl31|] ] 
[Pich and SC H ] 



M p = 762.5 + 2 MeV, 

M p = 754 ± 18 MeV, 

M p = 763.0 + 0.2 MeV, 

M p = 764.1 + 2.7^;° MeV, 



T p = 142 ±7 MeV, 
T p = 148 ±20 MeV, 
T p = 139.0 + 0.5 MeV, 

T p = 148.2+ I.9I59 MeV. 



(28) 



5. Conclusions 

The Pade approximants are important tools for the analysis of QCD amplitudes. They provide alternative deter- 
minations with competitive precision. 

At low energies, we have been able to determine the first derivatives of the 7T7t-VFF at the origin, this is, its 
quadratic charge radius ( r 2 (= 6a i) and the curvature Cy (= ai). For this, we used Euclidian data, i.e., from q 2 < 0. 

Likewise, we were able to employ the Pade approximants to extract information about the hadronic resonance 
poles (the p(770) mass and width). The experimental Minkowskian VFF and 7T7r-scattering data were analyzed by 
means of Pades centered between the first and second production thresholds. We obtained the determinations 
M p = 763.7 ±1.2 MeV and T p — 144 ± 3 MeV, with a competitive precision compared to other more elaborated and 
complex methods, in spite of the simplicity of the proposed procedure. 

This study shows the Pade approximants, once again, as a useful tool for the investigation of QCD phenomenology. 
They provide alternative determinations and, in spite of their simplicity, they have been proven as an efficient and 
systematic instrument. Nevertheless, when no theorem supports the convergence of the Pade sequence, the extracted 
parameters may have little to do with the physical ones, as we saw in ourLSM study of "Pade"-unitarizations in Sec. 3. 
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